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ABSTRACT 


Objective: The objective is to obtain the stresses due to strip loading in orthotropic plate lying over an irregular isotropic elastic 
medium. 

Methods: Anti-plane strain problem with perfect bonding boundary conditions following by Fourier Transformation on the equi¬ 
librium equation are used to obtain the solution. 

The deformation field due to shear line load at any point of the medium consisting of an orthotropic elastic layer lying over an 
irregular isotropic elastic medium is obtained. The anti-plane strain problem with the presence of rectangular irregularity is con¬ 
sidered. In order to study the effect of irregularity present in the medium and of anisotropy of the layer, we computed shearing 
stresses in both the media graphically. 
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INTRODUCTION 

It is well known that the upper part of the Earth is reeognized 
as having orthorhombie symmetry. Orthorhombie Symmetry 
is also expeeted to oeeur in sedimentary basins as a result of 
eombination of vertieal eraeks with a horizontal axis of sym¬ 
metry and periodie thin layer anisotropie with a vertieal sym¬ 
metry axis. When one of the planes of symmetry in an or¬ 
thorhombie symmetry is horizontal, the symmetry is termed 
as orthotropie symmetry and most symmetry systems in the 
Earth erust also have orthotropie orientations (Crampin^). 

The problem of deformation of a horizontally layered elastie 
material due to surfaee loads is of great interest in geoseiene- 
es and engineering. In material seienee engineering, the ap- 
plieations related to laminate eomposite material are inereas- 
ing. Many works related to Earth, sueh as fills or pavements 
eonsist of layered elastie medium. When the souree surfaee is 
very long, then a two-dimensional approximation simplifies 
the algebra and one ean easily obtain a elosed form analyti- 
eal solution. The deformation field around mining tremors 
and drilling into the erust of the Earth ean be analyzed by the 
deformation at any point of the media due to strip-loading. It 
also eontributes for theoretieal eonsideration of voleanie and 


seismie sourees as it aeeount for the deformation fields in the 
entire medium surrounding the souree region. It may also 
find applieation in various engineering problems regarding 
the deformation of layered isotropie and anisotropie elastie 
medium (Garg et aP, Singh et aP). 

The study of statie deformation with irregularity present 
in the elastie medium due to eontinental margin, mountain 
roots ete is very important to study. Chattopadhyay"^, Kar et 
aP, De Noyer^, MaE, Aeharya and Roy^ diseussed the prob¬ 
lems with irregular thiekness. Love^ provided the solution of 
statie deformation due to line souree in an isotropie elastie 
medium. Salim^® studied the effeet of reetangular irregularity 
on the statie deformation of initially stressed and unstressed 
isotropie elastie medium respeetively. The distribution of the 
stresses due to strip-loading in a regular monoelinie elastie 
medium had been studied by Madan et aP\ The effeet of 
rigidity and irregularity present in fluid-saturated porous 
anisotropie single layered and multilayered elastie media on 
the propagation of Love waves had been analyzed by Madan 
et aP^ and Kumar et aP^ respeetively. Reeently, Madan and 
Gabba^"^ studied two-dimensional deformation of an irregular 
orthotropie elastie medium due to normal line load. 
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In this paper, we have obtained the elosed-form expressions 
for the displaeement and shearing stresses in a horizontal or- 
thotropie elastie plate of an infinite lateral extent lying over 
an irregular isotropie base due to strip-loading. Numerieally, 
at different sizes of irregularity, we have studied the varia¬ 
tions of shearing stresses with horizontal distanee and it has 
been observed that the shearing stresses show signifieant 
variation with horizontal distanee at the different depth lev¬ 
els. 


with eo-ordinate planes as planes of elastie symmetry are 


(Jj + ^2^2 

2 ^^ 21^1 ^ ^^ 22^2 ^ ^^ 23^3 

^3 “ ^ 13^1 ^ 23^2 ^ 33^3 

^23 “ ^^44^23 
^13 “ ^^55^13 
^12 “ ^^6^12 


( 6 ) 


PROBLEM FORMULATION 

Consider a horizontal orthotropie elastie plate of thiekness H 
lying over an infinite isotropie elastie medium with x^-axis 
vertieally downwards. The origin of the Cartesian eoordinate 
system (x^, x^, X 3 ) is taken at the upper boundary of the plate. 
The orthotropie elastie plate oeeupies the region 0<x^<H 

and is deseribed as Medium I whereas the region x^> H is 
the isotropie elastie half spaee over whieh the plate is lying 
and is deseribed as Medium II. (Fig. I) 


where are normal strain eomponents and ^ 12 ?^ 23’^13 

are normal strain eomponents. The suffiees c.. {i,j = 1 , 2 , 3 , 4 , 5 , 6 ) 
are stiffnesses of an orthotropie elastie material. 

The strain - displaeement relations are given as 


1 

^ 12-2 


' dui 

.6x2 dxi) 


and ei 


dx 


etc. 


(7) 


In terms of displaeement eomponents, the equilibrium equa¬ 
tions ean be written from equations (3) - (7) as : 


Suppose a shear load P per unit area is aeting over the strip 

I X 2 |< /? of the surfaee x^ = 0 in the positive x^-direetion. The 
boundary eondition at the surfaee x^ = 0 is 

\-P\x,\<h 


The irregularity is assumed to be reetangular with length 2a 
and depth d. The equation of the reetangular irregularity may 
be represented as 


= f/(x2) = | 


J I X 2 |< a 

0 I X 2 |> a 


( 2 ) 


where £ = — « 1 is the perturbation faetor. 
2a 


6 +F 12 ) 


^2 6 x 16 x 2 


^66 


6x1 


(9) 


'"^-+(C44 + C23):""^ 


6 x 16 x 2 


6 x 16 x 2 


+ C55^+C44^+C33^=0 


( 10 ) 


Consider the field equation of an orthotropie material in anti 
- plane strain equilibrium state as: 


= W 2 = 0 , W 3 = W 3 (Xi,X 2 ); 


( 11 ) 


The non-zero stresses for an anti - plane strain equilibrium 
state are 


THEORY 

The equilibrium equations of Cartesian eoordinate system 
(x^, X 2 , X 3 ) for zero body foree are 


7'i2,2 7’i 3,3 — b 

(3) 

1 ^2,2 ^' 23 ,3 ~ ® 

(4) 

1 7'32,3 ^3,3 “ ® 

(5) 


_ _ ^ dU3 

""31 - ^55 ^ 


"^32 — Q- 


6 x 2 


( 12 ) 

(13) 


Equilibrium Equations for an orthotropie elastic medium due 
to anti - plane strain deformation are found to be 


6^U2 


m 


2 ^ ^3 


= 0 


(14) 


where m = 



where <j^,<J 2 ,<j^ are normal stresses and ^i 2 ?^i 3 ’^ 2 p^ 23 ’^ 3 p ^32 
are ealled shearing stresses. 

The stress-strain relations for- an orthotropie elastie medium 


At the interface (y, x = sf(y)), the boundary conditions are: 

1 . u ’3 =u‘-i . 

2 . - iefXy)z'32 = rsi - ief'iyWsz- 


(15) 
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By using the boundary condition (15), we find the deforma¬ 
tion field at any point of the orthotropic elastic plate corre¬ 
sponding to irregular contact between the orthotropic plate 
and the isotropic elastic half space due to strip-loading. 

Taking the Fourier transform of the equilibrium equation 
(14), we get 

^ 0 (16) 

dxl \Css y dxi Css 

The solution of the ordinary differential equation is 

u‘^ = ( 17 ) 

where and may be functions of k 
By using inverse Fourier transform, we have 

“3 = (18) 

Using equation (12), (13) and (18), the shear stresses are 


T 32 = ^[—im + €26 ^^^'‘kdk] 


2n L CO 


( 20 ) 


-/ ■ 1 u 

T 31 =- sm kh 


Therefore 


From (19) and (21), we obtain 


Cl-C2 = 


2P 


( sin kh\ 
k\k\ ) 


( 21 ) 


( 22 ) 


(23) 


■^32 = ——C^e 


(26) 


T(k' — msf 'Cy')^ _ 7 €26 + (fc' + sf XyS) ^2^ 

(28) 

where T = — and k' = 


Solving (23), (27) and (28), we get 


Tifelfel 


'[k'{V-e-‘ 


(fc'+£/'(y)K')e-^^"^l'^l^lyl) 


■2£m |fe|/(y))-e/ '(y)K'(H-e“2£m |fc|/(y)) J 

(k'+Ef'(y)V')e-^^ |fe|/(y)) 


Tik\k\ V" ' fe'(r-e-2^"'l''l/C>'))-e/'(y)r'(l+e-2^"'l''l/C>’))y 

2Psin kh ( k'{l+V) _ \ ^-e(m-l)\k\f(y) 


T^k\k\ 


{^k'{V-e-‘ 


■2£m |fe|/(y))_£^'(y)K'(l+e“2^'" 1^1/(y ))J 


where V = (T - 1)/(T + 1). 


(29) 

(30) 

(31) 


PROBLEM SOLUTION 

By applying Fourier Transformation technique on equation 
(2) we obtained 


/(/c) = -^sin(/ca) 


(32) 


Where = mCss = VQt^^ss- U- Using the boundary condition 
(1), we have 


Therefore, by using inverse transformation, we have 

f(y) = ^sin(/ca) dk = sign{a — x^) + sign{a + X 2 ) (33) 

where is the signum function. 

By substituting the values of constants and Q from 
equations (29), (30), (31) in the equations (18), (19), (20) 
for Medium I and in (24), (25), (26) for Medium II and also, 
substituting the value of f(y) for rectangular irregularity from 
equation (33), we obtain the following expressions for dis¬ 
placement and stresses. 


Ti V k\k\ 

The displacement in the isotropic elastic half space > H is 
uH = J_r°® (24) 

■J 7 tt *'— go ^ 


For Med. I 

r(l + yngrn\k\(2nE{signia-X2)+sign(a+X2)))'^ ^gm|fe|xi ^ 

y'g-m|fe|xi)|g-ifex2 


(34) 


(1 + V) tan- 

I 


IXi \ 

xl-h^J 


+ > V” tan 


The coefficient is to be determined from the boundary 
conditions. 

From equations (12), (13) and (17), we obtain 

(25) 


X2 + m^x\ 

2hm(2ns{sign(a - X 2 ) + sign(a + X 2 )) + Xi) 
xf + w?-{2ne{sign(^a - X 2 ) + sign(a + X 2 )) + Xi)^ - 
2hm(2n£(^sign(a - X 2 ) + signia + X 2 )) — Xi) 
xf + m^(2ne{sign(a - X 2 ) + sign{a + X 2 )) — Xj)^ — 


. Pm 
^32 = -t- 


(l + K)log 


log 


(h + X 2 )^ + m^xi 
(/i - X2)2 + m^xl 

{h + X2Y + m?( 2 nE{sign{a - X2) + sign{a + X2)) + Xi)^ 


Equations (15), (18), (19), (20), (24), (25) and (26) yield the 
relation 

C^e^rn\k\f(y) + |/c|/(y) ^ C2e-^\^\f^y^ (^7) 


+ r log 


(h - X2)2 + m^{2n£{sign{a - X2) + sign(a + X2)) + Xi) 
(/i + X2)^ + m?{2ne(sign{a — X2) + sign{a + X2)) - Xi)' 

{h — X2)2 + m^{2nE(sign{a — X2) + sign{a + X2)) - Xi)' 


(35) 


(36) 
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For Med. II 


P r sin kh ^ / 

^ ^|l:|(ta+l)4.isn(a-„)+«sn(a+„)) -Xl)g-ifeX2^/^ 


(37) 


= -^(1+1') 
7rri 


-I- 


2ft (( m + l)£(^sign(a — X2) + sign{a + X2)) — 
xf + (pn + l)e(signia - X2) + sign{a + X2)) - 
2 h (( 2 m(n + 1 ) + l)e(sign(a - X2) + sign(a + X2)) - x^ 
xf + (^( 2 m(n + 1 ) + l)E(sign(a - X2) + sign{a + X2)) - 


(38) 


Pg 


+ V) 


log 


(/i + X2)^ + (pn + l)e(signia - X2) + sign(a + X2)) - Xi^ 

(/i - X2)2 + (pn + l)e(signia - X2) + sign(a + X2)) - Xi^ 

^ (h + X2)^ + f( 2 m(n + 1 ) + l)8(sign(a - X2) + sign(a + X2)) - Xi) 

+ X log-^- h 

(/i - X2)2 + {p2m(ji + 1 ) + l)e{signia - X2) + sign(a + X2)) - Xi j 


(39) 


NUMERICAL RESULTS AND DISCUSSION 


In this section, we intend to examine the effect of irregularity 
on the stresses due to shear line load acting at any point of 
the orthotropic elastic layer lying over an irregular isotropic 
half space. For numerical computation, we use the values of 
elastic constants of Topaz (Orthotropic) for Medium I and 
the values of elastic constants of Glass (Isotropic) for Me¬ 
dium II given by Love^. 

Figures (2)-(4) and Figures (5)-(7) show the variation of 
shearing stresses and respectively, with horizon¬ 
tal distancefor different values of = 1,1.2,1.4,1.6 and for 
different depth levels = 0.5,1,1.5 . Figures (5)-(7) clearly 
show that for different values of, the difference between 
shearing stresses in magnitude significantly decreases as the 
depth increases. 

Figures (8)-(10) and Figures (11)-(13) show the variation of 
shearing tI\ and respectively with horizontal distance¬ 
for different values of a = 1,1.2,1.4,1.6. It has been found 
from the Figures (8)-(10) that for different values of a, the 
difference between shearing stresses in tI\ magnitude sig¬ 
nificantly increases as the depth increases. 


CONCLUSIONS 

The explicit expressions for the shearing stresses in an elas¬ 
tic medium consisting of orthotropic elastic layer lying over 
an irregular isotropic half space due to shear loading has 
been obtained. The results obtained are useful to study the 
static deformation around mining tremors and drilling into 
the crust of the Earth. The results are also useful to study the 


effect of irregularity present between the layer and the half¬ 
space. Graphically, it has been observed that the difference 
between the shearing stresses in magnitude in orthotropic 
elastic layer decreases as depth increases due to irregularity 
present. 

Further, it has also been observed that in isotropic semi- 
infinite half-space, the difference between the stresses in 
magnitude increases with the increase of depth. Thus, it has 
been concluded that the stress distribution in a layer with ir¬ 
regularity present at the interface is affected by not only the 
presence of irregularity but also by anisotropy of the elastic 
medium as a result of shear load acting over the strip of an 
orthotropic elastic medium. 


ACKNOWLEDGEMENT 

Authors acknowledge the immense help received from the 
scholars whose articles are cited and included in references 
of the manuscript. The authors are also grateful to authors/ 
editors/ publishers of all those articles, journals and books 
from where the literature for this article has been reviewed 
and discussed. The authors are also extremely thankful to the 
reviewers and editors for helping in the improvement of the 
paper. 


REFERENCES 

1. Crampin S. Geophysical Prospecting 1989; 37: 753-770. 

2. Garg NR, Madan DK, Sharma RK. Two-Dimensional Deforma¬ 
tion of an Orthotropic Elastic Medium due to Seismic Sources. 
Phys. Earth Planet. Inter., 1996; 94: 43-62. 

3. Singh K, Madan DK, Goel A, Garg NR. Two-Dimensional Static 
Deformation of Anisotropic Medium. Sadhana(India) 2005; 30: 
565-583. 

4. Chattopadhyay A, Chakraborty M, Pal AK. Effect of irregularity 
on the propagation of guided SH-waves. J. Meehan. Theor. Ap¬ 
plied 1983; 2: 215-225. 

5. Kar BK, Pal AK, Kalyani VK. Propagation of love waves in an 
isotropic dry sandy layer. Acta Geophysics 1986; 157:157-170. 

6. Noyer JD. The effect of variations in layer thickness of Love 
waves. Bull. Seism. Soc. Am. 1961; 51: 227-235. 

7. Mai AK. On the frequency equation of love waves due to abrupt 
thickening of crustal layer. Geofis. Pure Applied 1962; 52: 59- 
68 . 

8. Acharya DP, Roy I. Effect of surface stress and irregularity of 
the interface on the propagation of SH-waves in the magneto¬ 
elastic crustal layer based on a solid semi space. Sadhana 2009; 
34: 309-330. 

9. Love AEH. A Treatise on the Mathematical Theory of Elasticity. 
Dover Publication, New York 1944. 

10. Selim MM. Effect of Irregularity on Static deformation of Elas¬ 
tic Half Space. International Journal of Modem Physics 2008; 
22:2241-2253. 

11. Madan DK, Chugh S, Singh K. Stresses in an anisotropic Elastic 
Plate due to Strip-Loading. International Journal of Mechanics 
2011; 5: 57-62. 

12. Madan DK, Kumar R, Sikka JS. Love wave propagation in 


Int J Cur Res Rev | Vol 9 • Issue 4 • February 2017 















Madan et.al.: Stresses in an orthotropic elastic layer lying over an irregular isotropic elastic half-space 


an irregular fluid saturated porous anisotropie layer with rigid 
boundaries. Journal of Applied Seienee and Researeh 2014; 10: 
281-287. 

13. Kumar R, Madan DK, Sikka JS. Shear wave propagation in mul¬ 
tilayered medium ineluding an irregular fluid saturated porous 


stratum with rigid boundary. Advanees in Mathematieal Physies 
2014; 10: 1-9. 

14. Madan DK, Gabba A. 2-Dimensional Deformation of an Irregu¬ 
lar Orthotropie Elastie Medium. lOSR Journal of Mathematies 
2016; 12:101-113. 


Med. / 


H 


Med. If 





Orthotropie 


Isotropic 



Figure 4: Variation of the Shearing Stress in Med. I with 
the horizontal distance at x^= 1.5. 
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Figure 1: Section of the Model. 




Figure 5: Variation of the Shearing Stress ^2 in Med. I with 
the horizontal distance X 2 at x^= 0.5. 


Figure 2: Variation of the Shearing Stress in Med. I with 
the horizontal distance X 2 at x^= 1. 




Figure 6: Variation of the Shearing Stress ^2 in Med. I with 
the horizontal distance X 2 at x^= 1. 


Figure 3: Variation of the Shearing Stress in Med. I with 
the horizontal distance X 2 at x^= 1. 
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Figure 11: Variation of the Shearing Stress 


Figure 7: Variation of the Shearing Stress 7^2 in Med. I with the horizontal distance at x^= 0.5. 



Figure 8: Variation of the Shearing Stress tI\ in Med. II with the horizontal distance X 2 at x^- 1. 
the horizontal distance X 2 at x^= 0.5. 



Figure 9: Variation of the Shearing Stress in Med. II with 
the horizontal distance X 2 at x^= 1. 



Figure 10: Variation of the Shearing Stress in Med. II with 
the horizontal distance X 2 at x^= 0.5. 
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